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We study the quantum partition function of non-relativistic, ideal gas in a (non-cubical) box falling 
freely in arbitrary curved spacetime with centre 4- velocity Ubox- Using perturbed energy eigenvalues 
to evaluate the canonical partition function, we find that corrections to various thermodynamic 
quantities such as mean energy, entropy and specific heat include a very specific, sub-dominant term 



characterized by the dimensionless quantity, H 



i?ggA 



where _Rn 



BOX '^BOX 



and A — (She. 



This H-contribution does not depend on kinematic details of the system such as box dimensions and 
mass of particles, and in particular leads to Ss = (1/2)/3L'h (see text), a relation familiar from black 
hole thermodynamics. What is curious is that our result depends crucially on quantum mechanics 
since, in effect, the gas is allowed to "feel" the presence of the box through use of unperturbed wave 
function satisfying appropriate boundary conditions at the box walls. This is the feature which a 
classical analysis will completely miss. Our result might bear relevance in several contexts, such as 
analyses related to the generalized second law, a better understanding of horizon thermodynamics, 
and the "emergent gravity" viewpoint. 
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I. INTRODUCTION 



In the light of several intriguing connections between 
gravity and thermodynamics discovered over the past few 
years (see [1] for a recent review), it is perhaps useful to 
go back and study thermal properties of "normal" sys- 
tems (such as ideal gas) in presence of gravity [2]. The 
importance of this cannot be overemphasized, since it 
is possible that certain features of black hole thermody- 
namics are simply features of standard thermodynamic 
systems when curvature of spacetime is accounted for in 
the analysis of the latter [3]. But of course we know 
that black hole thermodynamics is inherently quantum 
mechanical in origin, and hence one may not learn any- 
thing "drastically" new from curvature corrections to 
standard systems using classical statistical mechanics; 
at best, classical analysis would yield terms represent- 
ing tidal forces (needed to hold the box together) etc. It 
is more interesting to ask whether a quantum mechani- 
cal calculation can yield any new information. It is this 
particular question which we address in this note, in the 
context of one of the simplest thermodynamic systems - a 
box of ideal gas. We consider such a box of ideal gas in an 
arbitrary curved spacetime, with its center freely falling 
along a geodesic Mbox, and compute corrections to the 
partition function due to spacetime curvature. And in- 
deed, we find that all thermodynamic quantities acquire a 
specific correction term (besides others) which is indepen- 
dent of system details such as box dimensions and mass 
of particles. This term is characterized by the dimen- 



sionless quantity S = R^A , where i^gg = RabU% 



^BOX'^BOX 



and A = I3hc. We highlight several features of this 
contribution, in particular the fact that S-b. = (1/2)/3J7e:, 
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a relation which is familiar from black hole thermody- 
namics. It would be worthwhile to investigate other pos- 
sible implications of these corrections from the viewpoint 
of thermodynamic aspects of gravity. 

We describe the set-up and relevant calculations in the 
next section; the perturbed free-particle energy eigenval- 
ues are given in Eqs. (6), and the final results are given 
in Eqs. (9). In the final section, we conclude with a 
few remarks on implications of the result. In order-of- 
magnitude arguments, we will use TZ to denote typical 
magnitude of curvature tensor components. 



II. THERMODYNAMICS OF BOX OF IDEAL 
GAS 

Consider a box of gas, whose center is falling freely in 
a curved spacetime with 4- velocity ixbox- The spacetime 
in a local neighbourhood of this trajectory can be con- 
structed using Fermi normal coordinates, in which the 
metric takes the form [4] 

m = - [(1 + a^y'f + Rd.s.y'y"] + ^iy') 
55. = -(2/3)i?o.,.2/V + o(y') 

5.. = S^, - (1/3) i?.p.a2/''2;" + 0{y^) (1) 

where a = Vu^^j^Mbox = in our case. In these coordi- 
nates, Mbox = 9q (i.e., the original trajectory is simply 
2/^ = 0), and we define a "box" as a confined region with 
flat "coordinate" faces, i.e., —{Lx/2) < x < -\-{Lx/2) 
and similarly for y and z. This box is filled with an ideal 
gas at temperature /3~^. The Hamiltonian for the con- 
stituent particles can be written in Fermi coordinates as 
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H = — Pq, which can be motivated as follows. The action 
for a single particle can be written in covariant form as 
J Pidx"^. In Fermi coordinates, the time coordinate as- 
signed to each point is the same as the proper time along 
the trajectory Mbox, say r. It is then natural to define 
3-velocities of the particles with respect to t. Hence, the 
action can be split as / (pg + p^dx'^/dr) dr. Comparing 
this with the usual form of the action, we can immedi- 
ately identify H = — pg. Note that this Hamiltonian is 

different from H = p • Mboxi where itfox is u^ox parallel 
transported to the location of the particle. This point is 
crucial for interpretation issues, and we will come back to 
it in a moment. Wc can write out the Hamiltonian explic- 
itly in terms of the metric coefficients using p-p = —m^; 
this gives 



H = 
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00 
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Let us now get back to the point concerning physical 
significance of this Hamiltonian. Of course, choice of 
a Hamiltonian depends on the time coordinate along 
which the Hamiltonian generates the evolution. By def- 
inition of Fermi coordinates, this time coordinate is r, 
and hence H = p • dr- In particular, since the metric 
is rjab all along the trajectory in Fermi coordinates, this 
Hamiltonian correctly describes the additional "gravity" 
field which would appear for particles located away from 
the trajectory, taking the zero of gravity potential at 



the trajectory. To illustrate the point, compare this 
with the case of a static box in Schwarzschild, with its 
center at some r ~ ro- After some simple manipula- 
tions, it is easy to see that the Hamiltonian as con- 
structed above becomes, in the non-relativistic limit, 
H — mc^ = /2m + [GM /r'^)my, which is the correct 
expression when the gravity potential at r = is chosen 
to be zero. In fact, when the spacetime has a timelike 
Killing vector, then the so called "energy at infinity" is 
what the above Hamiltonian gives; in this sense, the role 
of 

i^Box — dr above is the same as that of the Killing 
vector when available. These conceptual issues are also 
significant for what we mean by temperature T = 
What /3~^ represents is the temperature of a reservoir 
(say, a heating element) at the center of the box, which 
is our reference point. This /3 is therefore constant, unlike 
the Tolman temperature which is defined with respect to 
local energy of the particles. The above discussion is im- 
portant to correctly describe a thermal equilibrium state 
as a Gibbs state in a curved spacetime. Another way to 
emphasize the physical aspect of this discussion is that, 
since heat has weight, a quantity of heat flowing from one 
point to another point at a different gravity potential will 
be changed in content by its gravitational potential en- 
ergy, and hence the correct temperature assigned at a 
higher potential is lower than the temperature at a lower 
potential by the Tolman redshift factor. 



In the non-relativistic limit, the Hamiltonian reduces 



to 



^ " = 1^ + ^o.y'y'^) + ^^c^^g^o.y'y^ + \R',.,,y''y''cp, - l^^p^^y'v"^-^ (3) 
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Further, in the non-relativistic limit, terms which go as 
TZ-y^ X (p/mc)^ are second order of smallness, and there- 
fore can be ignored. Moreover, it follows from Einstein 
field equations that TZ ^ Gpm/ where pm is mass den- 
sity. Using all this, it is easy to see that, in the c — oo 
limit, we are finally left with 

H-mc^ = ^ + \^c'R-,,'^Xy'' (4) 

Incidentally, this expression matches with the non- 
relativistic limit of Dirac Hamiltonian obtained by Parker 



(sec Eq. (9.13) in [5]), which provides further separate 
support for the arguments leading to it. In what follows, 
we shall ignore the time dependence carried by curva- 
ture components. This is a reasonable assumption if time 
scale on which curvature changes, TZ/TZ, is much larger 
than typical time scale associated with the gas; we expect 
this to be the case at high enough temperature. 

The energy eigenvalues can be easily found using first 
order perturbation theory, with the unperturbed eigen- 
functions 
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corresponding to different parities. Here V — L^LyL^ is the box volume and ??i (l,oo). We assume that 
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the sides of the box are incommensurable, so that non- 
degenerate perturbation theory can be used. Then one 
obtains 

° 1 f 6 \ 



1=1,3 * 



(6) 



fundamental significance - indeed, it is this term which 
will turn out to have an interesting form. 

The canonical partition function turns out to be 



InZ 



In Zf + PNmc 
9'(0) 



g"(o) 



87r2 



+ E 



47r2 



The partition function for N particles, assuming Boltz- 
mann statistics, is Z = /Nl where z is the one par- 
ticle partition function. This can be calculated by ap- 
proximating intermediate sums as integrals, and assum- 
ing /V <C 1, where A = h/\/2TrmkT is the thermal de 
Broglie wavelength. First, define the functions 

n=l 

q{s) = sp{s) (7) 

The behaviour of q{s) can be deduced by properly ap- 
proximating sums with integrals. However, as we shall 
see, the final results can be stated simply in terms of q{s) 
and its derivatives at s = 0, so all that is required is that 
q{s) be analytic at s = 0. Also, in what follows, terms 
of the form "RX'^/L" have been ignored, since they are 
small compared to the terms retained by factors of X/L. 
The only sub-dominant term retained below is the one 
which is independent of box dimensions and mass m 
of constituent particles, and hence expected to be of some 



where \nZp = ln(y^A^'^^/A^!) is the flat space expres- 
sion. To put the final results in a neat form, it is conve- 
nient to introduce the following 

Definitions: 

A = /3hc which is a length scale independent of 
mass TO (unlike A), and is therefore more funda- 
mental 

^« = and S, = L,/X (> 1) 

ci = -(7"(0)/(27r) and = g'(0)/(27r); numerical 
values of these are not relevant, but can be shown 
to be ci = 1/2, C2 = 7r/12 

C4orr ^ U — Up and S'corr = S — Sp , whcrc 
Uf = 3N/2(3 and Sf = 3N/2 + N\n (eV/NX^) are 
standard flat space expressions 

It is now straightforward to use standard definitions 
U = —8/3 In Z and S = In Z + l3U to evaluate Ucon, S'corr 
and heat capacity at constant volume, Cy ~ —0^dfjU ~ 
3N/2 + CvcoTT- We obtain 



25corr/iV = +Cii?ggA2 - CsA^ ^ TZ.S, 

1=1,3 

l3Ucorr/N = +Cii?ggA2 - (1/24) CsA^ ^ TZA {36 - 2nc^^ 6^) 

i=l,3 

Cvcorr/N = -CiR^A^ + (3/4 )c2A2 (9) 

1=1,3 



which are correct to 0{5^ ). These are the final expres- 
sions wc wanted to derive, with several features worth 
highlighting: 

1. The Li and to independent part of the corrections, 
which we denote by S, depend only on S = R^^^ 

with = Rab Mbox ""box- 

2. Further, 5h and t/s satisfy 

5h = i/3C/H (10) 



where S~. = (A^ci/2)S and so on. This is Euler re- 
lation of homogeneity two, and is well known from 
black hole thermodynamics; in particular, black 
hole horizons have temperature, entropy and en- 
ergy which also satisfy S = {l/2)/3U. 

In fact, such a relation between S and U acquires 
importance in the emergent gravity viewpoint, since 
it can be interpreted as equipartition law for energy 
of microscopic degrees of freedom associated with 
spacctime horizons [6] . It is therefore quite intrigu- 
ing that the S correction terms, which appear due 
to quantum mechanics and are independent of sys- 
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tem details, have features in common with thermal 
features of spacetime horizons. 

The relevance of such Euler relation and area scal- 
ing of entropy for self-gravitating systems has also 
been emphasized in [7]. 

3. The S contribution to specific heat is negative if the 
strong-energy condition {Rqq > 0) holds. 

Incidentally, note that it is possible to incorporate finite 
size effects in the calculation rather simply by approxi- 

oo 

mating ^ exp(— an^) ~ (y^n/a — l)/2 instead of just 

n=l 

^jT^lal^ as is usually done; if one does this, the resultant 
corrections turn out to be of the same "form" as above 
with certain additional terms; in particular, the entropy 
correction gets an additional term of the form cq^/A^ 
where cg = N^? jV and A = 2{LxLy + LyL^ + LzLx) is 
the total area of the box. 

Finally, let us point out an important fact concern- 
ing the flat spacetime limit of the result. The curva- 
ture terms which appear in Fermi coordinates refer to the 
background spacetime, say S, which must be evaluated 
on the trajectory Mbox which is a geodesic of S. How- 
ever, since our system has finite energy, for consistency 
we must consider a geodesic in the spacetime compris- 
ing of background and the box contents. 5 4- where 
B is the perturbation to S caused by the box contents. 
Therefore, even when S is flat, there will always be some 
contribution from curvature produced by energy (rest -I- 
thermal) in the box. 



III. DISCUSSION AND IMPLICATIONS 

It is well known that thermal behaviour of systems 
with long-range interactions exhibit several peculiar fea- 
tures, and their statistical analysis requires considerable 
amount of care [8] . Some of these features such as neg- 
ative specific heat and unequal "local" temperature in 
thermal equilibrium became more widely known since 
the discovery of thermodynamic aspects of black holes. 
However, these same aspects, which have their origin in 
quantum field theory in non-trivial coordinate systems, 
also indicate that one must take a closer look at standard 
thermodynamics in an external gravitational field. For 
example, temperature and entropy associated with local 
acceleration horizons are observer dependent, and one 
must therefore try to find a natural way to incorporate 
these inevitable features in conventional thermodynam- 
ics. These and related issues form the main motivation 
for the analysis in this note. 

Let us first summarize the main steps of the calcu- 
lation. We defined the Hamiltonian appropriately using 
Fermi coordinates, reduced it to the non-relativistic limit, 
and then moved over to quantum mechanics by using first 
order perturbation theory. We must emphasize that the 



effect of the box is rather subtle, and is not just a pure 
finite size effect, since the 5 terms are independent of Li 
and m. 

One context in which the result might be relevant is the 
so called generalized second law (GSL), which was for- 
mulated by Bekenstein [9] based on the fact that horizons 
have entropy (proportional to area). Quantum mechan- 
ics played a crucial role in Bckenstein's analysis, since 
classically the minimum increase in entropy of horizon, 
when a particle (or a system) falls across it, is zero. In 
our case, entropy of the freely falling box of ideal gas ac- 
quires the additional term proportional to S = /JqqA^ 
(besides other terms depending on details of the sys- 
tem). This term appeared because we used quantum 
mechanical wave functions to evaluate the perturbed en- 
ergy eigenvalues and hence the partition function. How- 
ever, note that the curvature correction in Eq. (6) docs 
not have h in it - the only non-trivial aspect at this level 
is the dependence of energy eigenvalues on 1/nf. Never- 
theless, quantum nature of the result is evident in the S 
corrections which are 0{H^) since A = 0{h). It would be 
interesting to investigate how (if at all) does this affect 
the analyses related to GSL. 

It is, of course, important to generalize the result pre- 
sented here to different systems, and also to the rela- 
tivistic case. However, the fact that the relevant term 
depends only on the length scale A = I3hc suggests that 
the result will survive in the relativistic limit (since it 
does not depend on mass) . ^ A rigorous analysis of other 
systems (such as harmonic oscillator) confined in a box 
is complicated by the fact that even unperturbed energy 
eigenvalues and eigcnfunctions arc not known in analytic 
form. One hopes to arrive at a better understanding of 
deeper physical implications of the result provided such 
technical issues could be overcome. Incidentally, note 
that the time scale At = A/c = H/kT is the time un- 
certainty associated with thermal fiuctuations of energy 
kT. There is a hint here of interplay between thermal and 
quantum fluctuations, but again, more work is needed to 
understand the underlying physics better. 

I thank T Padmanabhan for discussion. This research 
is funded by NSERC of Canada, and Atlantic Association 
for Research in the Mathematical Sciences (AARMS). 



Appendix A: Heat Capacity and Average Pressure 

Another important observation can be made regarding 
average pressure of the gas. Pressure on the i—th face 
can be defined as (3Pi = {Li/V)d\nZ/dLi since V/Li = 
A±i is the transverse area of the face. If one defines 



In fact, a close look at the results in [5] indicates that relativistic 
corrections produce qualitative changes only in Z ^ energy levels 
of Hydrogen atom, which strengthens this belief further. 
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the correction to average pressure by P = (1/3) Pi = 
■Fcorr+^F, where /3Pf = N/V is the flat space expression, 
then one can easily show that 

As mentioned above, S-b. = ^/^f/s- Hence, what the 
above expression says is that, the average pressure is con- 
tributed by those parts of ^corr sjid Ucon which do not 
satisfy S = (1/2) f3U. ^ Using the above result, one can 



write down the correction to equation of state of the gas: 
^ = 1 - (t) E ^'^^ (12c2 + 27Td,) (A2) 

^ ^ j=l,3 

Note that the final pressure does depend on geometric 
features of the box. In fact, it might be possible to in- 
terpret the second term in the round brackets on RHS as 
additional pressure due to tidal forces acting on the box. 
Quantum mechanics adds to this the C2 piece (note that 
A^^, - O(fti) while A^5f - 0(ft°)). 
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